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I. INTRODUCTION 



Statistics plays a central role in physics, and marks one of the main parting points 
between the classical and quantum domains. In conventional phenomena, taking place in 
3 + 1-dimensions, the invariance under the interchange of identical particles forces the wave 
function to be either symmetric or antisymmetric. The spin-statistics theorem in relativistic 
quantum field theory then ties this symmetry to the spin, and all fundamental excitations are 
classified as either fermions or bosons. Through the spin the fermionic or bosonic character 
can manifest itself even in one- or few-particle systems. However, it is in the thermodynamics 
of a system that a symmetric or antisymmetric wave function has a dramatic effect. 

In 2 + 1-dimensional phenomena new possibilities open up, and excitations of generalized 
statistics interpolating between fermions and bosons may occur. The role of these anyons |TJ 
in the quantum Hall effect is fairly well established 0, and they may be relevant also in 
the superconductivity of some materials at high temperature 0. Of course, in these and 
other possible phenomenological applications in two-dimensional condensed matter systems, 
anyons are just low energy descriptions of the effective behaviour of particles and interactions 
which, although restricted to two-dimensions, exist in three-dimensional space and conform 
to the conventional boson-fermion classification. 

However, the same reasons which make anyons possible at low energies, pose a challenging 
question. It is certainly of great theoretical interest to establish the extent to which anyonic 
excitations can exist as fundamental objects in a truly 2 + 1-dimensional world, just as 
ordinary fermions or bosons, or if on the contrary, they are always a low energy illusion. 
Mathematically, the possibility of generalized statistics has its origin in the topology of 
the configuration space of a many particle system in 2 + 1-dimensions. Under the double 
interchange of two identical particles, x\ — x 2 — > x 2 — X\ — > X\ — x 2 , the observability of the 
probability density requires that the wave function transforms as ip — > e l2e ip. In spacetime 
dimensions greater than three this double interchange is just the identity operation, so we 
are forced to 9 = (bosons) or 9 = 7T (fermions). In 2 + 1-dimensions however, the winding 
of the trajectories followed by the particles as they are interchanged may be non trivial, in 
which case the operation is not connected to the identity and 9 remains unrestricted. 

This, however, is only a mathematical possibility. Its realization in a quantum theory 
depends on non trivial issues regarding the short distance behavior of the theory. Indeed, the 
very notion of the linking of two trajectories assumes some kind of hard core repulsive inter- 
action between the particles (accounting for some generalized exclusion principle), and this 
may or may not be consistent with a formulation from first principles. In a non relativistic 
treatment, the transformation law of the anyonic wave function under particle interchange 
provides a well defined meaning for the concept of generalized statistics. In practice this 
is implemented by minimally coupling the matter fields to an Abelian Chern-Simons gauge 
field. However, the non existence of an ideal anyon gas, analogous to the Bose and the 
Fermi gases, represents a major obstacle in the understanding of anyonic thermodynamics. 
A major tool in this respect has been the perturbative expansion in the statistical angle 9, 
which when done at the bosonic end exhibits the singular nature of the statistical interaction 
at short distances Q. Although at the computational level this can be regularized with a 
repulsive 5-function contact potential, these short distance difficulties are beyond the scope 
of the non relativistic theory. 



2 



The physics of anyons at short distances must be formulated and discussed in the frame- 
work of relativistic quantum field theory. This has been done mostly in the canonical 
formalism |J, where the main difficulty is again the non existence of free asymptotic many 
particle states. Besides posing technical difficulties which are still controversial ||, this ob- 
scures somewhat the concept of generalized statistics in this framework. To a large extent, 
it has acquired in the literature a rather algebraic meaning, in terms of the phase factor 
which enters generalized canonical commutation relations, again interpolating between the 
conventional commutator and the anticommutator for bosons and fermions. But those alge- 
braic relations are rather formal, so long as there is no direct relation with free asymptotic 
states, and they do not address the relevant short distant issues, so they may or may not 
turn out to be consistent with a proper regularization procedure. 

Since one expects that the statistics obeyed by the fundamental excitations will have di- 
rect consequences in the thermodynamic properties, a possible way around these difficulties 
is to use the functional integral formulation to study the theory at finite temperature, thus 
avoiding such difficult issues as the lack of free asymptotic states or the exact kinematical 
meaning of generalized statistics in the relativistic theory. Here we adopt that point of view, 
and report our first results for some thermodynamic properties of massive fermions in three 
space-time dimensions coupled to an Abelian Maxwell-Chern-Simons field. In particular, we 
evaluate the specific heat at finite temperature and density for different values of the statis- 
tical angle 9. The calculation is done to leading order in the conventional loop expansion, 
complemented with a low momentum approximation which allows for further simplifications. 

We should admit from the start a limitation of our approach. At least within the loop 
expansion, the theory is not properly defined unless the Maxwell term is included in the pure 
gauge action. (Note, however, that a Maxwell term is generated anyhow by the quantum 
fluctuations of the matter fields, even if absent at the classical level). As it happens, that 
term is also expected to conceal the effect of the generalized statistic induced by the Chern- 
Simons term, at least at short distances ||. Thus, we do not expect to see in our results a 
very notorious shift, say, from fermionic to bosonic properties as we vary the statistical angle 
9. To complicate matters further, in 2 + 1-dimensions there is no Bose-Einstein condensation, 
so the specific heat does not exhibit a dramatic difference between fermions and bosons, as 
we are used to in 3 + 1-dimensions. Nevertheless, our results do show that the dependence of 
the specific heat in the statistical angle is consistent with a smooth transit from a fermionic 
to a bosonic behavior. 



II. THE THEORY AT ZERO CHEMICAL POTENTIAL 

We consider massive fermions in three Euclidean dimensions coupled to an Abelian 
Maxwell-Chern-Simons gauge field with a Lagrangian 

1 e 2 

£ = + m)i/j + -Fl v - i—e^xF^Ax + e</>4^, (1) 

where e 2 has dimensions of mass and the statistical angle 9 is a dimensionless parameter. 
Our choice of 7- matrices is ji = i<Ji, with <7j the usual Pauli matrices. 
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In order to establish conventions, describe our treatment and review previous results, we 
shall evaluate first the partition function at finite temperature and zero chemical potential. 
In the following section we introduce a non-vanishing chemical potential to discuss the 
thermodynamics of the system at finite fermion density. 

We carry out the finite temperature calculations as usual, compactifying the (Euclidean) 
time variable into the range 0<r</5=l/T(in our units, h — c — k — 1). Then, 
the functional integral defining the partition function should be computed using periodic 
(antiperiodic) boundary conditions (in time) for bosons (fermions). The partition function 
is then defined as 

Z = M(J3) J VyDtPDAp exp (S[A,}) , (2) 

with 

S[A^\ = dr J d 2 xC = J d 3 x C. (3) 

The relevance of the normalization factor N(f3) will be discussed below. Integrating out the 
fermions one has 

Z = H{{3) det ( _)(# + m) J VA^ exp [- (S[AJ + S q [AJ)] , (4) 

where the subindex (— ) means that the determinant has to be evaluated using antiperiodic 
boundary conditions, and 

S , K ]=-l D det H (l + ^_) (5) 

is the contribution of the fermionic quantum fluctuations to the effective action. To one 
loop order this is 

SM,} = ttr( 7 J—] 2 . (6) 



2 + 



m 



When computing the trace, the antiperiodic boundary conditions are implemented in mo- 
mentum space replacing integrals over p° by sums over the discrete Matsubara frequencies 



(2n + l)?r 



A = (7) 

The one-loop result for S g in Eq. (|) can be evaluated in closed form 0. However, 
to simplify the numerical analysis, here we shall work to leading order in an expansion in 
powers of the momenta. Then, following Ref. || one finds 

S g [A,} = h(P) d 3 x [^-Fl - —e^xF^Axj , (8) 

where the function h(P) is given by 
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■tanh(m/3/2). 



(9) 



Then, in this approximation the partition function is given by a Gaussian integral over the 
gauge field, times the partition function of free massive fermions: 

Z = M{(3) det H (i# + m) f VA^ exp (-S eff [A„]) , 



where the effective action S. 



S + S q is given by 



SefAM = \ j/x { (l + \ HP)) + 2mhC9)) e^ A i^A A j . (10) 

From the coefficient of the Maxwell term, one sees that, as usual, the loop expansion re- 
quires a = e 2 /47rm <C 1. Similarly, from the coefficient of the Chern-Simons term one would 
conclude that in this case the loop expansion also requires 9 <C 2n. Yet, from the expe- 
rience gained with three-dimensional fermions at zero temperature PJIDfl, one could hope 
that this restriction over 9 may not really apply. This, however, is controversial at finite 
temperature [TT| . 





As discussed in ||12|| , at finite temperature one has to take into account the contribution 
from the Faddeev- Popov determinant arising from gauge fixing, even if the theory is Abelian. 
We shall choose the Lorentz gauge (d^A^ = 0) for which the Fadeev-Popov determinant is 
given by 



A 



FP 



det {+) (-d 2 



(11) 



where now the (+) subindex indicates that the determinant has to be evaluated using peri- 
odic boundary conditions |T2j . Evaluating the contribution to the determinants coming from 
the space-time indices and keeping track of all temperature dependent factors, one obtains 



Z = det R (0 + m) x det^ 1/2 (-d 2 + M 2 ) 



where M 2 is defined by 



M 



e 2 /9 + 4mh((3) 
2 + (4/3)h(/3) 



(12) 



(13) 



and the subindex R denotes that one must take the finite part of the determinants: temper- 
ature dependent divergences coming from these determinants and from the normalization 
factor N{f3) cancel each other [fT^ . Also, following |i~3| , in Eq. ( |T2|) we have omitted an 
overall divergent factor which does not contribute to the thermodynamic properties of the 
system. 

It is interesting to note that the partition function (|T2| ) corresponds to an effective theory 
with a free massive fermion and a free boson with a mass M which is related to the Chern- 
Simons topological mass. This mass acquires a temperature dependence which arises from 
the fermion determinant. It should be stressed that all dependence on the statistical angle 
9 occurs through that mass. 



5 



Determinants are easily evaluated. Standard finite temperature methods |12| lead to 

d 2 p 



In detptiw + m ) 
lndetp(-«9 2 + M 2 ) 



V 
2V 



■ln(l +e" 



(2tt) 2 

|_0_ ln(1 _ e -/3v^; 

+2(3VAn R , 



(14) 



where AQr is the finite contribution to the bosonic determinant coming from the "zero 
point energy": 



An, 



d 2 p 



p 2 + M 2 



R 



12tt 



-M 3 . 



The thermodynamic potential is defined as 

1 



n 



(3V 



InZ, 



so that in the present case we have 



Q = — 



(15) 



(16) 



(17) 



(3 J (2vr)2 ^ ' (3 J (2tt) 

It should be stressed that we have not included in the thermodynamic potential the zero 



point energy contribution AQr, which is eliminated by an appropriate normal ordering |L4 
Now, using C v = -{3 2 d 2 (f3n)/d(3 2 , we get 



dx 



x 



+ 



4tt(3 2 Jpm cosh 2 (x/2) 8ir(3 2 Jpu sinh 2 (ic/2) 



dx ■ 



x 



1 dM 



P 2 M 2 



7T d(3 \exp(/3M) 

2 / 



+ 



1 d 2 M 
2^~d(3 2 ' 



! Mln(l - exp( 



+ 



1 dM 



-[3 3 M + j3 2 ln(l - exp(-pM)) - 



(3 3 M 



exp(-/3M) 



(18) 



Except for being at zero chemical potential, this result was the goal of our calculation. It 
is not usual in field theory to look at purely thermodynamic functions like C v . But, as this 
calculation has illustrated, for our purposes is important since it gives direct information 
about the quantum excitations of the system, without need of facing difficult issues related 
to the asymptotic states. In the next section we extend this calculation to incorporate a 
finite fermionic density. 



III. THE THEORY AT NON-ZERO CHEMICAL POTENTIAL 

To describe the system at finite fermion density we introduce a chemical potential /i, so 
the Lagrangian in Eq. ([I]) becomes 
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1 e 2 

C = + m- ilon)tl> + -F 2 U - i—e^\F^ v A x + eipfi^j. 
The contribution of the fermionic quantum excitations to the effective action is now 



lndet ( _) ( 1 + - 



up + m — ijon 



(19) 



(20) 



Comparing this with Eq. (§) we see that the presence of the chemical potential fi amounts 
to a shift of the Matsubara frequencies p° n in Eq. (^)) to 



(2n + l)vr . 

Pn = n ^ 







(21) 



Then, repeating the previous calculation with these shifted frequencies, a tedious but 
straightforward calculation leads to the effective action 



SefAAp] = i jT d 3 x I (l + ~ g(/3, //)) F 2 U -i^+2m g(/3, e^F^A^ 



where now 



16nm 



'P 



P, 



tanh — (m — fi) + tanh — (m + /i) 



(22) 



(23) 



This result holds, again, to leading order in a low momentum expansion. Moreover, we have 
neglected in this expression Lorentz non-invariant terms which are of order /3exp(— /3m). 
These terms are negligible for f3m sufficiently big, which is anyway assumed by the low 
momentum approximation. 

Following the same steps as in Sec. [H], for the thermodynamic potential we obtain 



n = -- 



i 



i 



(3 J (2ny 

d 2 p 



ln(l + e" 



d 2 p 



(3 J (2tt) 2 



ln(l 



(3 J (2tt) s 



ln(l + e -/3(Vp 2 + m2 -M) ) 



where now 



e 2 /# + 4mg(/3,/z) 
2 + (4/3)g(/?,/i) 



(24) 



(25) 



For the specific heat we now find 
1 



ln(3 



§(m+ti) 

1 dM M ( 



(x — Pfi)x 2 
cosh 2 (z/2) 



, (x + du)x 2 
dx — — — + 



0{m-ii) 



cosh 2 (x/2) 



dx 



x 



sinh 2 (x/2) 



(3 2 M 2 



(/') 



1 «9 2 M, 



7T d(3 \exp(PM w )-l 
dM M \ 2 ' 



2?r <9/3 2 



/3 2 M (At) ln(l-exp(-/3M (M) )) 



1 



2tt V 9/3 



ft + P 2 ln(l - exp(-/3M M )) - 



P 3 M ifl) 



exp(-(3M {fl) ) 



(26) 
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Similarly, for the fermion number density p = 



—dn/d/j,, we get 



d 2 p 

i 



1 _|_ e P{y/p 2 +m 2 +Li) i _|_ e P{^/^+r 



2tx(3 



-Mi 



dM< 



hi 



2sinh^M M 



(27) 



This is a constraint equation which must be inverted to determine the chemical potential 
p in terms of a given density and temperature. It should be noticed that the first two 
terms in this expression correspond to the free fermion (relativistic) gas, with both fermions 
and anti-fermions contributing to the density. Similarly, in Eq. (|2q) , the first three terms 
correspond to the contributions to C v from free fermions and anti-fermions of mass m and 
chemical potential ±/z, and from free bosons of mass Mi^ without a chemical potential. 
The numerical analysis of these equations is the subject of the next section. 



IV. NUMERICAL RESULTS 

In this section we present our numerical results. For given values of the statistical angle 9 
and the fermionic density p, we invert Eq. (|27j) numerically (with a precission of one part in 
10 5 ), and then substitute the chemical potential into Eq. (p6|) to obtain the specific heat C v 
as a function of the temperature at different values of p and 9. Since here both e 2 and m have 
dimesions of mass, we find it useful to parametrize the theory in terms of the dimensionless 
variables a = e 2 /4irm, p/m 2 , (3m and 9. 

We have performed numerical computations covering the following domain of parameters: 

a E (HT 2 , KT 6 ) 

4 e ( 10 "M°) 

#G(l(T 4 ,i) (28) 



and in the temperature range 



(3m G (1, 10 2 ) (29) 



There are no qualitative differences for the results over the whole domain of parameters. 
Typical curves for the specific heat as a function of the temperature are presented in Figs. [I] 
and | for 9 = 1 and 9 = 1(T 4 , respectively, at a = 9.73xlCT 7 , and p/m 2 = 1.78 xl(T 8 . These 
are values relevant for high T c superconductors |T5 |, and correspond to m — 7.5 x 10 10 cm -1 , 



thrice the value of the bare electron mass, and a fine structure constant in 2 + 1 space-time 
dimensions obtained from that in 3 + 1 space-time dimensions dividing by the interplanar 
spacing. 

Although the curves in Figs. [I] and |2] seem almost identical, the specific heat in Fig. [I 
(9 = 1) is slightly larger than that in Fig. |2] (9 = 10~ 4 ). This is a characteristic feature, 
which we find to hold in the entire parameter range, although the actual difference depends 
on the values of a and p. We think it is a very interesting result, since it is consistent 
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with the (unimpressive) difference between the free bosonic and fermionic specific heats in 
2 + 1-dimensions. This is illustrated in more detail in Figs. § and ||. In Fig. [3] we present the 
difference of the specific heats C v (bosonic) — C v (fermionic) for the free relativistic (charged) 
boson/fermion gases, both at p/m? = 0.01. This is to be compared with Fig. f|, where 
we present the difference C v {9 = 1) — C v (9 = 10~ 4 ), both at a = 0.01 and p/m 2 = 0.01. 
Although this is not a proof of a shift from fermionic towards bosonic statistics (whatever 
that means in this relativistic system), it is certainly consistent with that intuitive picture: 
since we start from a fermionic system at 9 = 0, one expects a transition towards a bosonic 
behavior as 9 is increased, provided the nonrelativistic ideas about Chern- Simons-generalized 
statistics hold in this case. 

In connection with the curves shown in Figs. |l|-§, it is worthwhile noticing that their 
apparent violation of the classical Dulong-Petit law is not so. For a 2 + 1-dimensional non- 
relativistic classical ideal gas one finds C v = p. In our case, however, p stands for the Abelian 
charge density, and therefore is the density of particles minus the density of anti particles, 
as shown by the first terms in Eq. fl27|) in the case of fermions. Hence, one can indeed reach 
large values of C v without a corresponding increment of p. 

We have studied also the fermion density as a function of the chemical potencial. We 
obtained a smooth function which does not exhibit the step structure found in |TB| where it 
is interpreted as a signal of a critical behavior. Note, however, that that reference considers 
a high temperature range (opposite from ours) and starts from a pure Chern-Simons system, 
without a Maxwell term. 



V. CONCLUSIONS 

We have analyzed at finite temperature and density a system of massive fermions in three 
space-time dimensions coupled to an Abelian Maxwell-Chern-Simons field. In particular, we 
evaluated the specific heat for different values of the statistical angle 9. This was done to 
leading order in the loop expansion and in the low momentum approximation. In doing 
so, we have shown that it is possible to obtain physical information of the system without 
having to give an explicit answer to such difficult issues as the lack of free asymptotic states 
or the exact kinematical meaning of generalized statistic in relativistic quantum field theory. 

Our main result, the specific heat as a function of temperature and density (Eqs. 



27)), was analysed in a wide range of parameters within the low-temperature regime. This 
limitation to low temperatures arises in part from the low momentum expansion, which we 
here adopted for simplicity but can be avoided 0. However, to move to higher temperatures 
one will also have to face the contribution of Lorentz non-invariant terms which here we have 
discarded. We hope to come back to this point in a future work. 

As we ourselves had expected, the specific heat has only a mild dependence on the 
statistical angle 9. However, this dependence is consistent with what one expects from the 
nonrelativistic picture of generalized statistics, thus suggesting that those ideas may still be 
applicable, to some extent, in the relativistic domain. 

As we pointed out earlier, a more marked dependence on the statistical angle 9 may be 
expected in a pure Chern-Simons theory without Maxwell term. That, however, will possibly 
require a non-perturbative treatment. In this respect it may be of interest to consider other 
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thermodynamic functions which may be more sensible to the long-distance properties, and 
therefore more immune to the presence of the Maxwell term. 
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FIGURES 

FIG. 1. Specific heat as a function of temperature at a = 9.73 x 1CT 7 , p/m? = 1.78 x 10~ 8 and 
= 1. 

FIG. 2. Specific heat as a function of temperature at a = 9.73 x 1CT 7 , p/m 2 = 1.78 x 10~ 8 and 
6 = 10" 4 . 

FIG. 3. Difference of the specific heats C„(bosonic) and (^(fermionic) for the free relativistic 
(charged) boson/fermion gases, both at p/m 2 = 0.01. 

FIG. 4. Difference of the specific heats C v (d = 1) and C v (6 = 10~ 4 ), both at a = 0.01 and 
p/m 2 = 0.01. 
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